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In this work, motivated by the energy bound suggested by Cohen et ai, we propose the modified 
holographic dark energy (MHDE) modeL Choosing the IR cut-off L — Rcc and considering the 
parameterizations = 2 — Xa, = 2 — 3Aa^/ (l + 3a^) and — 2 — Aa^/ (/3 + a^) , we derive all the 
physical quantities of the non-saturated MHDE model analytically. We find that the non-saturated 
MHDE models with the parameterizations = 2 — Aa and n'^ — 2 — 3Aa^/ (l + 3a^) are single- 
parameter models in practice. Also, we consider the cosmological constraints on the non-saturated 
MHDE, and find that it is well consistent with the observational data. 
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I. INTRODUCTION 



Recently, the holographic dark energy (HDE) has been considered as an interesting candidate of dark 
energy, which has been studied extensively in the literature. It was proposed from the well-known 
holographic principle [2] in the string theory. For a quantum gravity system, the local quantum field 
cannot contain too many degrees of freedom, otherwise the formation of black hole is inevitable and then 
the quantum field theory breaks down. In the black hole thermodynamics 0, |3| , there is a maximum 
entropy in a box of size L, namely the Bekenstein- Hawking entropy bound Sbh, which scales as the 
area of the box ~ L^, rather than the volume ~ L'^. To avoid the breakdown of the local quantum 
field theory, Cohen et al. 5] proposed a more restrictive bound, namely the energy bound. If pA is the 
quantum zero-point energy density caused by a short distance cut-off, the total energy in a box of size L 
cannot exceed the mass of a black hole of the same size namely pi^L^ < nipL, where nip = {8nG)~^^^ 
is the reduced Planck mass. The largest IR cut-off L is the one saturating the inequality. Thus, 

PA - Sn^mlL-^, (1) 

where the numerical constant is introduced for convenience. If we choose L as the size of the 
universe, for instance the Hubble horizon H^^, the resulting p\ is comparable to the observational 
density of dark energy 0]. However, Hsu [Ti] pointed out that in this case the resulting equation-of- 
state parameter (EoS) is equal to zero, which cannot accelerate the expansion of the universe. The other 
possibility [1] is to choose L as the particle horizon 

Rn^aj -^a I — , (2) 



a 



where H = d/a is the Hubble parameter; a = (l-|-z)~^ is the scale factor of the universe; z is the redshift; 
we have set ao = 1; the subscript "0" indicates the present value of the corresponding quantity; a dot 
denotes the derivative with respect to cosmic time t. However, it is easy to find that the EoS is always 
larger than —1/3 and also cannot accelerate the expansion of the universe in this case with L = Rh [9|. 
To get an accelerating universe, Li proposed in Q to choose L as the future event horizon 



oo 



dt dd 



a 



Rh = a I - = a I —. (3) 



In this case, the EoS of the holographic dark energy can be less than —1/3 [9]. In fact, the version of HDE 
with L ~ Rh is the one considered extensively in the literature. For a comprehensive list of references 
concerning HDE, we refer to e.g. , [III, m, m, [11] and references therein. 

In the literature, many authors extensively considered the HDE whose energy density was taken to be 
the one given in Eq. ([1]). Here, we instead propose to consider a modified version of HDE. Notice that the 
energy bound proposed in requires that the total energy in a box of size L cannot exceed the mass of 
a black hole of the same size. Therefore, it is reasonable to argue that when HDE coexists with matter, 
the total energy should satisfy the energy bound, namely, ptotL'^ ^ m^L. Correspondingly, the largest IR 
cut-off L is the one saturating the inequality. Thus, 

Plot = 3n^mpL"^, (4) 

where ptot = Pm+PA is the total energy density; pm = Pmoa~^ is the energy density of matter. Therefore, 
the energy density of the modified HDE is given by p\ = ptot — Pm = i'n?rn^L~'^ — Pmoo"'^ instead. Note 
that we do not consider the radiation in this work, since almost all the cosmological observations probe 
the cosmic history after the matter-radiation equality. Using the Friedmann equation Sm^iJ^ — ptot and 
Eq. Q, it is easy to find that H^LF' = n^. Note that we consider a flat universe throughout. Since n is 
given in the form of throughout, we only consider the positive n in this work. So, we have 

HL ^ n. (5) 

This very simple equation is our starting point for the modified HDE (hereafter MHDE). Naturally, the 
next step is to choose an appropriate L. 
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II. CHOOSING AN APPROPRIATE IR CUT-OFF 



The most natural choice for L is L = H~^. In this case, one should have n = 1 from Eq. ([5]). So, the 
requirement of the holographic principle, namely Eq. (U, is nothing but the familiar Friedmann equation. 
There is no new input. Therefore, the choice L = is trivial. 

The second choice for L is the particle horizon Rh in Eq. ([2|). For L — Rh, we have 

da a Ha a \ n J 

which can be regarded as a differential equation for L with respect to a. From Eq. one can find that 
L cx a^+" ^ and H = n/L oc a~^~" \ Noting that in the matter-dominated epoch H cx a~^^^, MHDE 
with L = Rfj can describe the matter-dominated epoch when n ^ 2. This is a good news. However, we 
should consider the second thing. As is well known, 

^ H + H^ = H^ + Ha— . (7) 

da 



a 



Substituting H oc a^^^" into Eq. ([7]), we have d/a — —n^^H^ < 0. Therefore, the universe cannot be 
accelerated, and hence the choice L = Rh is unsuitable. 

The third choice for L is the future event horizon R^ in Eq. For L = Rh, we have 

^ = ^-J- = L(l-i]. (8) 
da a Ha a \ nj 

So, it is easy to see that L oc a^~" ^ and H = n/L oc a^^+" \ Unfortunately, it cannot describe the 
matter-dominated epoch in which H oc a~^^^. On the other hand, substituting H oc 0"^+" ^ into Eq. ([7]), 
we have a/a = n~^H'^ > 0. Thus, the universe is accelerating, but cannot undergo a decelerated epoch. 
So, the choice L = Rh is also unsuitable. 

There is another choice. In [l^, Gao et al. suggested that L can be proportional to the Ricci scalar 
curvature radius. The resulting HDE was called Ricci dark energy in the literature. See e.g. [3, HI, [H, [13 
for works on the so-called Ricci dark energy. In [l^ [3, [iBj EB, Ii3| , there is no physical motivation to this 
proposal for L in fact. Recently, in 18|, Cai et al. found that the Jeans length Rcc which is determined 
by Rqq = H + 2H^ gives the causal connection scale of perturbations in the flat universe. Since the 
Ricci scalar is also proportional to H + 2H^ in the flat universe, Cai et al. in fact found the physical 
motivation for the Ricci dark energy. Here, we take L = Rcc as the fourth choice for MHDE. In this 
case, from Eq. ([5]), we obtain 

H'^ = n^L-^ =n^ {h + 2H'^^ . (9) 

It can be recast as a differential equation for H with respect to a, namely 

dH / 1 \ H 

-d^-[^-V^- (10) 

So, wc find that H oc a" When — 2, MHDE with L = Rcc can describe the matter-dominated 
epoch in which H oc a~^^^. On the other hand, substituting H (x a^ ~^ into Eq. ([7|), we find that 

'^^(^^i)h^. (11) 

Thus, the universe is decelerating when n > 1, whereas the universe is accelerating when n < 1. So, the 
choice L — Rcc gives the possibility for a cosmologically feasible MHDE model. 

However, in the saturated MHDE model, n is constant. For a given n, it is impossible to describe 
both the decelerated and accelerated phases, even in the case with L = Rcc- To build a cosmologically 
feasible MHDE model, we should go further. 
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III. NON-SATURATED MHDE 



In the derivation of HDE mentioned in Sec. HI the holographic bound was chosen to be saturated. 
However, this is not necessary. In [l9| . Guberina et al. proposed the so-called non-saturated HDE, in 
which the parameter n = n{t) is a function of cosmic time t, rather than a constant. Similarly, here we 
consider the non-saturated MHDE with L = Rcc, in which n = n{a) is a function of scale factor a. To 
solve Eq. pop . we should specify the explicit form of n{a). Similar to the well-known EoS parameterization 
w — wo + Wa{l~ a) — wi—W2a considered extensively in the literature we parameterize — a — Xa, 
where a and A are constants. In fact, one can regard this type of parameterization as a linearized 
expansion with respect to scale factor a. As shown in Sec. HIl when = 2, we find that MHDE with 
L = Rcc can describe the matter-dominated epoch in which H cx a~^^^. So, we get a = 2 by requiring 
— > 2 when a —> 0. From now on, we consider the non-saturated MHDE in which L = Rcc a-nd 

2 - Xa. (12) 

Noting that > 0, we should require A < 2 at least to describe the cosmic history up to now, namely 
< z < oo. Substituting Eq. into Eq. (fTU]) . we can recast it as 

(ia(^2 — Aa)a ^ ^ 

The solution reads 

H^Hq{2- X)^'^ (2 - Xay^'^ , (14) 

where Hq — H{z = 0) is the Hubble constant. Obviously, H oc a^^^^ when a 0, and hence this model 
can describe the matter-dominated epoch. For convenience, we rewrite Eq. (|14p to 

E' = {2-X){l + zf(^2-^^ \ (15) 

where E = H/Hq. Substituting Eq. into Eq. ((T]), we can obtain the deceleration parameter 

a 1 / Aa \ 

Obviously, the universe is decelerating when Aa < 1, whereas the universe is accelerating when Aa > 1. 
The transition occurs when Aa = 1. Therefore, the transition redshift is given by 

zt = X-l. (17) 

From the Friedmann equation irn^H^ = ptot = Pm + PA, we have pA = SrripH^ — pmoo,^^ , or equivalently 



The fractional energy density Qi = Pi/ (SnipH^^, where i = m and A. In particular, we find that 



2- A V 1 + ^, 

Requiring ilm — > 1.0 (namely matter dominated) when z oo, we have 

^rnO - 1 - ^ • (20) 

Unlike the other versions of HDE considered in the literature, r2„iO is not an independent parameter in 
the non-saturated MHDE model. This is an important result. Using Eq. ([^n]). we finally obtain 

^™ = i-2(iT^' (^1) 
i^A^i-a,.^^^. (22) 
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Interestingly, we find that CIa = at the transition redshift zt given in Eq. ()17|) which is determined 
by g = 0. Provided that A is of order unity, the cosmological coincidence problem could be alleviated. 
On the other hand, it is easy to find the effective EoS as (see e.g. [2ll. ISTj) 

Wcs = = -1 - -— = - 2g- 1 = -— — , (23) 

ptot 3 3 (2 - Aa) 

Noting that Wcff — QmWm + ^awa and Wm = 0, the EoS of HDE reads 

where f^A is given in Eq. (|22p . If A is given, all the physical quantities can be derived accordingly, thanks 
to the important result Eq. ((20|) . In fact, to our knowledge, the non-saturated MHDE model considered 
in this work is the fourth single-parameter cosmological model besides the well-known flat ACDM model, 
the flat DGP braneworld model [22J and the new agegraphic dark energy (NADE) model [23.] . 



IV. COSMOLOGICAL CONSTRAINTS ON THE NON-SATURATED MHDE 

Here, we consider the cosmological constraints on the non-saturated MHDE. In [U,!!^, Supernova 
Cosmology Project (SCP) collaboration released their new dataset of type la supernovae (SNIa), which 
was called Union compilation. The Union compilation contains 414 SNIa and reduces to 307 SNIa after 
selection cuts (24| . 

We perform a analysis to obtain the constraints on the sin gle p arameter A of the non-saturated 
MHDE model. The data points of the 307 Union SNIa compiled in [24| are given in terms of the distance 
modulus fiobs{zi). On the other hand, the theoretical distance modulus is defined as 

Mtft('Zi) = 51ogio£'L(2i) + Mo , (25) 
where no = 42.38 — Slog^Q h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 



The from the 307 Union SNIa are given by 



- ^) ' (27) 

where a is the corresponding la error. The parameter /iq is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over /ig. However, there is an alternative 
way. Following [1^ 113] , the minimization with respect to po can be made by expanding the xfi of Eq- l|27|) 
with respect to /j,o as 



2 



where 



xl{\)^A~2^ioB + ^,lC, (28) 



A(X) = V [^°''^(^'') ~ f^th{zt;tio = 0, A)]^ 
i ''la.s (^0 



^/x^x ^ Mofcs(g») - f^th{zi;p.Q = 0, A) C = V 
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Eq. ([28|) has a minimum for /uq — BjC at 

X^(A) = i(A)-^. (29) 



Since = obviously, we can instead minimize xi which is independent of /ip. Note that the 

above summations are over the 307 Union SNIa compiled in |24| . It is worth noting that the corresponding 
h can be determined by /io = BjC for the best-fit parameter. 

At first, we consider the cosmological constraints on the non-saturated MHDE, by using the 307 Union 
SNIa dataset only. In the top-left panel of Fig. [1] we present the likelihood £ oc e"-^ 1"^ versus A. We find 
that the best-fit model parameter is A = I.502too3o (with \a uncertainty) or A = 1.502^0 063 (with 
2a uncertainty), while Xmm = 311.272. The corresponding h = 0.707. We plot ii{z) for the best-fit 
parameter A in the top-right panel of Fig. [T] We also present the observational H{z) data from [11] 
(see also [2ll . [29l l30l |) for comparison. It is easy to see that H{z) of the non-saturated MHDE is well 
consistent with the observational i?(z) data. In addition, 57^, S^a, <?, fi'cff and wa as functions of redshift 
z for the best-fit A are also presented in Fig. [TJ Clearly, and fi™ are comparable since the near past, 
whereas fimo = 1 ^ A/2 = 0.249 and f^Ao — A/2 = 0.751. We see that q changed from g > to g < at 
Zt — \ — \ = 0.502. On the other hand, wa crossed the so-called phantom divide w = —1. 

In the literature, the shift parameter i? from the cosmic microwave background (CMB) anisotropy, 
and the distance parameter A of the measurement of the baryon acoustic oscillation (BAO) peak in the 
distribution of SDSS luminous red galaxies [3l|, are also used extensively in obtaining the cosmological 
constraints. The shift parameter R is defined by [H, 



where the redshift of recombination = 1090 which has been updated in the 5-year data of Wilkinson 
Microwave Anisotropy Probe (WMAP5) [131 ■ The shift parameter R relates the angular diameter distance 
to the last scattering surface, the comoving size of the sound horizon at z* and the ang ular scale of the 
first acoustic peak in the CMB power spectrum of the temperature fluctuations [H, Issj ]. The value of R 
has been updated to 1.710 ± 0.019 in WMAP5 [11]. The distance parameter A is given by 



1 Y" dz "1^/^ 



Zb Jq E{z) 



(31) 



where Zb = 0.35. In [H], the value of A has been determined to be 0.469 (ns/0.98)-°-3^ ± 0.017. Here the 
scalar spectral index Us is taken to be 0.960 which has been updated in WMAP5 [3^ . 

Now, we perform a joint x^ analysis to obtain the cosmological constraints on the single parameter A of 
the non-saturated MHDE model, by using the combined data of the 307 Union SNIa, the shift parameter 
R of CMB and the distance parameter A of BAO. Here, the total x^ is given by 

=XI+ XCMB + X%AO > (32) 

where x^ is given in Eq. (HH), Xca/b ^ {R- Robsf and x%ao ^ {A - AobsY I cr\- The best-fit model 
parameter can be determined by minimizing the total x^- In the top- left panel of Fig. [51 we present 
the likelihood C oc 1"^ versus A. We find that the best-fit model parameter is A = 1.461^q q29 (with 
Icr uncertainty) or A = 1.46lj;Q'Q6i (with 2a uncertainty), while Xmm = 320.858. The corresponding 
h = 0.699. We also plot B{z) for the best-fit A in the top-right panel of Fig. [21 Again, it is easy to see 
that ll(z) of the non-saturated MHDE is well consistent with the observational ll{z) data. In addition, 
rim, 5, Wcff and wa as functions of redshift z for the best-fit A are also presented in Fig.[2l Clearly, 
J7a and Sim are comparable since the near past, whereas rimo = 1 — A/2 — 0.269 and SIao = A/2 = 0.731. 
We see that q changed from g>Otog<Oatzt = A— 1 = 0.461. On the other hand, wa crossed the 
so-called phantom divide w — ~\. 

As shown above, the non-saturated MHDE with the simplest parameterization = 2 — Aa is well 
consistent with the observational data in fact. This is fairly impressive, since all the physical quantities 
of this model are given by the simple and analytical expressions. 
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FIG. 1: The top-left panel is the likelihood C ot e ^ ' versus A. The other panels are //, Q,m, Ha, q, Wes and 
wa as functions of redshift z for the best-fit A. These results axe obtained by using the SNIa data only. 
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FIG. 2: The same as in Fig. [T] except that these results are obtained by using the combined data of the 307 
Union SNIa, the shift parameter R of CMB and the distance parameter A of BAO. 
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V. OTHER PARAMETERIZATIONS FOR THE NON-SATURATED MHDE 



Although the simplest parameterization v? — 2~ \a given in Eq. has been shown that it is fairly 
successful to fit the cosmological observations, it is only valid for a < 2/A which is required by n? > 0. 
So, we cannot use it to describe the future evolution of the universe. In fact, it is non-trivial to find a 
suitable parameterization of n which is well-defined for the whole < a < oo. Naively, one can consider 
the simplest parameterization -n? ^ 2 ~ \a/ {(3 + a). However, we find that it cannot be consistent with 
the combined observational data of the 307 Union SNIa, the shift parameter R of CMB and the distance 



parameter A of BAO. Here, we consider another parameterization n 



), where a, (3 



and A are constants. Again, we obtain a = 2 by requiring — s- 2 when a 0, in order to describe the 
matter-dominated epoch in which H oc a^^^"^ . From now on, we consider the non-saturated AfflDE in 
which L = i?cc and 



= 2 



Aa2 



/3- 



(33) 



We require /? > to avoid singularity. Noting that < a^/ (/3 + oF) < 1 for < a < oo, we also require 
A < 2 to ensure > 0. Substituting Eq. ([33| into Eq. (fTO|) . we find that the solution is given by 



H = Hoa 



-3/2 



2/3 -f (2 - A) a2 



(34) 



Obviously, H oc a when a 
Then, we have 



2/3 + 2- A 

0, and hence this model can describe the matter-dominated epoch. 



On the other hand, 



= {l + zY 



2/3-f(2-A)(l + z)' 
2/3 + 2 - A 



(35) 



mO 



2/3 + 2 - A 



Again, requiring fi,, 



1.0 (namely matter dominated) when z 

X 



2/3 + (2 - A) (1 + z) 

oo, it is easy to find that 



^2/3 + 2- A, 

Therefore, Qmo is not an independent parameter. And then, we have 



2/3 



2(2-A) 



whereas JIa ~ I ~ Q,-, 



_2/3+(2-A)(l + z) 
Substituting Eq. into Eq. ([7]), we obtain the deceleration parameter 



am 



2/3 + (2 - A) a2 



(36) 



(37) 



(38) 



(39) 



It is easy to see that the universe is decelerating when (A — 1) < /3, whereas the universe is accelerating 
when (A — l)a^ > (3. The transition occurs when (A — 1) = (3. Therefore, A > 1 is necessary to 
accelerate the universe. The transition redshift is given by 



/A- 1 



(40) 
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The effective EoS is given by (see e.g. [2l|, [s? 



Ptot 
Ptot 



2 H 1 



1) 



3 [2/3 + (2 - A) a2 



whereas the EoS of HDE reads 



WA 



(41) 



(42) 



in which f2m is given in Eq. 

Similar to Sec. IIV( we perform a joint analysis to obtain the cosmological constraints on the model 
parameters A and /3, by using the combined data of the 307 Union SNIa, the shift parameter R of CMB 
and the distance parameter A of BAO. By minimizing the total given in Eq. (15^ . we find that the 
best-fit parameters are A = 1.9095 and /3 = 0.3455 (the corresponding h — 0.7024), while Xmin — 311.368. 
In Fig. [31 we present the 68% and 95% confidence level contours in the X — (3 parameter space. We also 
plot H, Urn, ^A, q, WeS and WA as functions of redshift z for the best-fit A and (3 in Fig. [H Obviously, 
the MHDE model with the parameterization given in Eq. ((55)) is well consistent with the observational 
data in fact. 
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FIG. 3: The 68% and 95% confidence level contours in the \ — (3 parameter space for the parameterization given 
in Eq. (|33p . The best-fit parameters are also indicated by a solid point. 



Note that there are two parameters in the parameterization given in Eq. (j33p . From the viewpoint 
of model building, it is better to find a single-parameter model. Inspired by the above results, we can 
consider another parameterization 



3Aa2 



(43) 



in which there is only one model parameter A. In fact, it is just the reduced version of Eq. (|33p with 
/3 = 1/3. So, all the physical quantities of the MHDE model with this parameterization can be obtained 
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FIG. 4: The H, f^A, Q, Weff and wa as functions of redshift z for the best-fit A and /3 for the parameterization 
given in Eq. (I33|l . These results are obtained by using the combined data of the 307 Union SNIa, the shift 
parameter R of CMB and the distance parameter A of BAO. 



by simply setting (i — 1/3 in Eqs. — fl^ . Fitting to the combined data of the 307 Union SNla, the 
shift parameter R of CMB and the distance parameter A of BAO, we find that the best-fit parameter 
reads A = 1.8843^^:^^3^ (with Icr error) or A = 1.8843;^ |^^g^ (with 2cr error), while = 311.405. The 
corresponding h = 0.7016. In Fig. [5l we present the corresponding Ukelihood C oc 1"^ versus A, as 
well as the i?, fim, f^A, 9, Wcff and wa as functions of redshift z for the best-fit A. Obviously, the MHDE 
model with the parameterization given in Eq. ()43|) is also well consistent with the observational data. 

Comparing to the parameterization (|12p . besides the advantage of being valid for the whole < a < oo, 
the parameterizations ((55)1 and (|15| have smaller Xmm when we fit them to the combined data of the 307 
Union SNIa, the shift parameter R of CMB and the distance parameter A of BAO. 




FIG. 5: The same as in Fig. (5] except for the parameterization given in Eq. (|43|l . 
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VI. CONCLUSION AND DISCUSSION 

In this work, motivated by the energy bound suggested by Cohen et al., we propose the modified 
holographic dark energy (MHDE) model. Choosing the IR cut-ofF L = Rqc and considering the param- 
eterizations ((HJ, (|33p and we derive all the physical quantities of the non-saturated MHDE model 
analytically. We find that the non-saturated MHDE models with parameterizations (jl2p and (j43|) are 
single-parameter models in practice. Also, we consider the cosmological constraints on the non-saturated 
MHDE, and find that it is well consistent with the observational data in fact. Given the simplicity and 
analyticity, the non-saturated MHDE model is very interesting and deserves further investigation. 

After all, there are some remarks on the parameterizations of n(a). Firstly, our situation is similar to the 
familiar parameterization for the EoS of dark energy considered extensively in the literature, such as the 
well-known w{a) = wo + Wa{l — a) Q and w{z) = wo+wiz (see e.g. [13, [H, H^l ) ■ Until we can completely 
understand the mysterious dark energy in the far future, the method of parameterization is useful to 
shed some light on the nature of dark energy. Secondly, as mentioned in Sec. IHI) our parameterization 

= 2 — Xa given in Eq. (|12p is motivated by the EoS parameterization w{a) = wq +Wa{l ~ a). However, 
as recently argued by Shafieloo et al. [4l|, the EoS parametrization w{a) = wq + Wa{l ~ a) cannot be 
trusted in the redshift interval < z < 1.4, and it fails for a ^ 1. Accordingly, one might doubt our 
parameterization — 2 — Xa. As shown in Sec. IIV( the parameterization — 2 ~ Xa works well for 
< a < 1, or, < 2; < cxD. However, as mentioned in the beginning of Sec. |Vl it is only valid for a < 2/X 
which is required by > 0. Thus, we proposed other parameterizations for in Eqs. p3p and (|43p 
which work well for the whole < a < cxo. Finally, one might argue that to be holographic (i.e., the energy 
density is proportional to the horizon area), n{a) should be nearly constant. Therefore, n{a) should not 
change so quickly. We admit that our parameterizations UH), ([55]) and (|43p do not fulfill this point. On 
the other hand, they can be used as working parameterizations since they are well consistent with the 
observational data in fact. Nevertheless, it is still worthwhile to seek a more satisfactory parameterization 
for n(a) which changes not so quickly and can be consistent with the observational data. We leave this 
to the future work. This tells us that there are still many improvements to do for the MHDE model. 
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